The electronic structure of the Na^CoC^ surface 
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The idea that surface effects may play an important role in suppressing e' g Fermi surface pockets 
on Na^CoC^ (0.333 < x < 0.75) has been frequently proposed to explain the discrepancy between 
LDA calculations (performed on the bulk compound) which find e g ' hole pockets present and ARPES 
experiments, which do not observe the hole pockets. Since ARPES is a surface sensitive technique 
it is important to investigate the effects that surface formation will have on the electronic structure 
of Nai/3Co02 in order to more accurately compare theory and experiment. We have calculated 
the band structure and Fermi surface of cleaved Nai/3Co02 and determined that the surface non- 
trivially affects the fermiology in comparison to the bulk. Additionally, we examine the likelihood of 
possible hydroxyl contamination and surface termination. Our results show that a combination of 
surface formation and contamination effects could resolve the ongoing controversy between ARPES 
experiments and theory. 



One longstanding controversial issue surrounding 
Na x Co02 is the discrepancy between between experi- 
mental measurements and theoretical predictions [l[ con- 
cerning the small hole pockets along the T— K direction. 
Angle resolved photoemission spectra (ARPES) exper- 
iments have clearly shown the existence of a large ai 3 
cylindrical hole on the Fermi surface which is consistent 
with LDA band structure calculations @, H, 0, 0| . How- 
ever, the small hole pockets which are also present in 
band structure calculations due to the partially filled e g ' 
band are absent. Whether this should be interpreted as 
a failure of LDA to account for some crucial correlation 
effects or as a failure of ARPES to detect the bulk struc- 
ture is a subject of much discussion. 

Dynamical mean field theory (DMFT) calculations 
that incorporate correlation effects on a local level and 
might be expected to bring theory and experiment into 
closer agreement, instead find that these only serve to 



the volume of the e s ' hole pockets 
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ertheless, since long range correlations such as ferromag- 
netic spin fluctuations are not treated by this method, 
one cannot rule out alterations to the Fermi surface 
caused by such correlations. 

Surface effects, as a broadly defined category, have also 
been proposed as a possible reason for the absence of 
the pockets in ARPES measurements. By probing the 
Na 2p core-level intensity as a function of photoelectron 
emission angle, it was clearly shown that strong differ- 
ences between the surface and bulk exist and that pho- 
toemission experiments contain a large, if not dominant, 
amount of surface character In particular, oxygen va- 
cancies are common on oxide surfaces allowing hydroxyl 
terminations to occur through a combination of oxygen 
vacancies and a water background in UHV or simply hy- 
drogen contaminants adsorbing on the surface [1 01 ] - In- 
deed, many previous studies of oxide surfaces have found 
that it is nearly impossible to remove water and hydrogen 
contaminants in a UHV environment 
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It is therefore desirable to understand, from a theo- 
retical point of view, the differences between the elec- 



tronic structure of a cleaved surface layer and that of the 
bulk material and to investigate the role of common de- 
fects such as unintended adsorbates on the fermiology of 
NaCo02- In this Letter, we use density functional theory 
(DFT) to examine the band structure and Fermi surfaces 
of a cleaved Na^OoC^ surface. We show that these dif- 
fer non-trivially on a surface as compared to the bulk 
and are sensitive to many uncontrollable factors, such as 
Na ordering, surface contaminants, and defects. Intrigu- 
ingly, we find that the surface is extremely susceptible 
to electron donor adsorbates and that these strongly de- 
press the e g ' hole pockets well beneath the Fermi energy, 
very similar to what has been observed in ARPES. 

Calculational Methods All calculations have been per- 
formed using a plane wave DFT code, the Vienna Ab- 
initio Simulation Package (VASP) [12j. The projected 
augmented wave method (PAW) [13( and the generalized 
gradient approximation [14] to the exchange correlation 
potential were used. The surface and bulk Na^OoC^ 
structures were modeled using a six stoichiometric layer 
slab (a 15A vacuum is added between successive slabs for 
surface calculations), so that both BZs are the same. For 
most calculations, we used a V3 x -\/3 hexagonal cell with 
three Co atoms per layer. The two top and bottom lay- 
ers were relaxed while the middle two were held in their 
bulk relaxed positions. Ground state geometries were 
fully optimized using a Monkhorst Pack k-point mesh of 
6x4x2. For total energy density of states (DOS) cal- 
culations, the Brillioun zone (BZ) integrations were per- 
formed on a 8x8x2 mesh. For some calculations, an even 
larger 2v / 3xv3 cell with six Co atoms per layer was re- 
quired. Na can occupy either the sites atop Co (Nal) 
or at the center of an equilateral triangle formed by Co 
atoms at the vertices (Na2), or a combination of the two. 
We have found that the energetics associated with Na2 
ordering are most favorable in the V3 x y3 cell, so we 
present those results here. This agrees well with exper- 
imental results from Roger et. al [la ] who observe from 
neutron diffraction experiments that when x=l/3, Na 
prefers the Na2 position, whereas at high doping, mul- 



2 



tivacancy clusters form with partial occupation of Nal 
sites. 

In order for the stoichiometry of the bulk compound 
to be maintained and for the surface Co atoms to have 
the same valency as those in the bulk, the number of 
Na atoms on each of the two surfaces should be half 
the number between bulk layers. For the x=l/3 com- 
pound, this requires one Na for every six Co and can- 
not be achieved with only 3 Co atoms per layer. We 
have therefore investigated an " over-stoichiometric" com- 
pound that contains one Na atom on each surface and 
an " under-stoichiometric" compound that has no Na on 
either surface. The stoichiometric compound was also 
investigated by further lowering the symmetry (doubling 
the cell size). Although we discuss the energetics of this 
calculation, the extensive backfolding of the BZ makes 
it impossible to visually interpret the Fermi surfaces. 
We have investigated the stability of the under- or over- 
stoichiometric surfaces by calculating the phase separa- 
tion energy difference El / 3 -+ El / 3 + > _g where Ei/ S _, 
-E/1/3+ and -E1/3 represent the energies associated with 
the under-stoichiometric, overstoichiometric, and stoi- 
chiometric systems, respectively. We find this inequal- 
ity to be true and therefore it is highly likely that the 
stoichiometric surface will not phase separate into areas 
that are either under- (Na depletion) or over- (Na excess) 
stoichiometric. 

We have adjusted the electron count in calculations of 
both under- and over-stoichiometric surfaces by shifting 
the Fermi energy appropriately to best mimic the most 
likely situation of a stoichiometric surface. As we dis- 
cuss, this approximation is not perfectly representative 
because the Na atoms in some cases have a non-negligible 
effect on the surface Fermiology. It does, however, allow 
us to "bracket" expected electronic structure modifica- 
tions between over- and under- stoichiometric results. 

Surface Calculations The bulk band structure and 
Fermi surfaces for Nai/3Co02 in a six-layer \/3 x \/3 unit 
cell is shown in Figs. [Iji and b. There are six sets of a\ g 
and eg Fermi surfaces (degeneracies prevent them from 
being individually distinguishable) , corresponding to the 
six different layers. Because of BZ backfolding the two 
types of Fermi level crossings occur along the same high 
symmetry line, T— M, in FigQJi and the e g '-derived small 
hole pockets are located inside the larger a\ g hexagonal 
pockets in FigQ}). For the bulk compound, locating the 
Na atoms at the Nal vs. the Na2 position has an imper- 
ceptible effect on the Fermi surface or band structure. In 
Fig. Hb (right panel) we show the Fermi surface for the 
over-stoichiometric surface compound. We present the 
over-stoichiometric case because although the effect of 
Na is exaggerated it more closely mimics the stoichiomet- 
ric case than neglecting it entirely. Moreover, the under- 
stoichiometric case also shows a noticeable decrease in the 
size of the e g hole pockets, though smaller than the over- 
stoichiometric case. The stoichiometric case, if it were 




FIG. 1: (a) Band structure of bulk and (b)(left panel) 
Fermi surface of bulk and (right panel) overstoichiometric 
Na!/ 3 Co02. The light colored solid bands show e 9 ' charac- 
ter. The oi 3 bands create the large, dark (black) cylindrical 
pockets. The elliptical hole pockets are shown in light (green) 



possible to calculate it in this unit cell, would lie some- 
where between the two. In the over-stoichiometric case 
one of the e g ' pockets is clearly shrunken and another has 
been suppressed entirely. These can be clearly identified, 
using atom-based projection, as the surface-originated 
pockets. The corresponding surface ai g bands, identified 
in the same way, expand proportionately. Thus, any tech- 
nique that probes the surface electronic structure cannot 
be expected to detect Fermi surfaces in agreement with 
bulk calculations and, additionally, surface sensitive re- 
sults should not be extrapolated to represent the bulk 
electronic structure. While the e g hole pockets are not 
fully eliminated by the existence of the surface, they are 
rather dramatically shrunken. The origin of this effect is 
in the noticeably shortened Na-0 distance (it decreases 
by 0.32 A at the surface). The Na+ Coulomb field low- 
ers the energy of the O electron bands and, by increasing 
the energy separation between O and Co states, decreases 
Co-0 bonding. Although we do not show it here, a pro- 
jection of the O character of the Co-derived t2 9 bands 
shows that Co-0 hybridization is stronger in the e g sym- 
metry parts of the band complex relative to the &i g parts 
and particularly strong at and near the part of the dis- 
persion that gives rise to the small e g fermi surfaces. 
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Thus, an energy lowering of the O states disproportion- 
ately affects these areas. 

In contrast to the bulk, the surface Fermi pockets (both 
small and large) are slightly sensitive to the Na position. 
The Fermi surfaces corresponding to Na in the Nal posi- 
tion show slightly larger e g hole pockets than when Na 
is in the Na2 position (shown in Fig. [TJd). The position 
of Na varies the size of the hole pockets due to electro- 
static interaction between Na + and the Co a\ g orbitals. 
The a\ g orbitals are shifted downward compared to the 
e s ' orbitals by positive Na + ions. This effect is stronger 
when Na is in the Nal position versus the Na2 position 
because the ai g orbital is elongated along the z-direction. 
The smaller shift of a\ g downward when Na is in the Na2 
position results in a smaller hole pocket 

The under- and over-stoichiometric compounds are, as 
noted previously, adjusted so that their electron count 
is equal to the stoichiometric one. Although the Fermi 
surfaces are mainly insensitive to this approximation, 
there is one notable difference that occurs due to the lack 
(under-stoichiometric) or surplus (over-stoichiometric) of 
Na on the surface itself. Na, being a free-electron metal, 
features a parabolic, nearly free-electron band as a free 
atom or when it serves as an intercalant in a compound. 
It was recently pointed that the position of this band 
in intercalated graphites is crucial for superconductivity 
[13]: in some compounds it is so low in energy that it 
crosses the Fermi level. It was observed that that the 
position of this intercalant band is quite sensitive to the 
interplanar distance [l^|. The interplanar distance de- 
termines how squeezed the Na electrons are between the 
host layers. Naturally, if the electrons are unbound on 
one side (surface), the intercalant band drops substan- 
tially (for CaCg this drop was found to be 0.5 eV [ljj]). 
In the bulk, the Na band is high enough to be neglected 
when considering the Co-0 bands, even at high Na con- 
centrations (it sinks upon adding Na, for electrostatic 
reasons). On the surface, however, it is 4.21 eV lower 
than in the bulk for the overstoichiometric compound. 
Fig. (2k- shows the overstoichiometric bandstructure in 
which the Na band drops down below the unoccupied 
e g complex and overlaps with the unoccupied a± g bands 
near the T point at approximately 250 meV above Ej?. 
For the stoichiometric compound (not shown), the bot- 
tom of the surface Na band is 670 meV higher than for 
the over-stoichiometric compound. Although the sur- 
face Na band does not cross E^ for the low Na content 
(2:= 1/3) considered here, it does dip substantially be- 
neath it for :e=2/3, and therefore changes the fermiology 
significantly. This may have relevance for the abrupt dis- 
continuity in hole count compared to Luttinger's theorem 
noticed in systematic ARPES experiments. 

The adsorption energy of a Na atom on the undcr- 
stoichiometric/stoichiometric surface is 3.92 eV/2.55 eV. 
These results indicate that it is highly favorable for Na 
to remain on the surface rather than desorb and that 



the cleaved surface of Na^OoC^ is highly susceptible to 
contamination by electron donor adsorbates. Although 
it seems unlikely that excess Na could be found to at- 
tach to the surface under UHV conditions, the presence 
of H2O molecules or H atoms is much more conceivable. 
Even assuming that a perfectly clean oxide sample was 
placed inside a UHV chamber, a H2O background and 
H contaminants are almost impossible to get rid of de- 
spite annealing at high temperatures and high vacuum 
pressures. In fact, it has been observed for some oxide 
surfaces that even after annealing to temps around 900K, 
H still appears to be stuck to the surface of the oxide 
[20l | . We have calculated that a single H atom will ad- 
sorb to the surface of the over-stoichiometric compound 
with an energy of 2.84 eV. For the stoichiometric/under- 
stoichiometric compounds this increases to 3.50 eV /4.12 
eV due to a reduction in Na-H repulsion. 

Surprisingly, the effect of adsorption of an H atom on 
the band structure and fermiology of Na^OoC^ is con- 
siderably more dramatic than simple electron donation 
would suggest. As can be seen by comparing the high- 
lighted e g surface bands of Fig. [5] a to the bulk e' g bands 
of Fig. [IJl, the presence of a H contaminant severely de- 
presses the entire e g ' band complex. The suppression is 
strongest along T-M where the bands formerly crossed 
Eip without the contaminant, but the downshift in en- 
ergy is obvious throughout the bandstructure. The ap- 
proximately -110 meV shift is congruent with what is 
seen in ARPES measurements [H, Hf, though no conclu- 
sive reports of surface contamination have yet been pub- 
lished. Fig. [2k illustrates the effect of H on the over- 
stoichiometric surface while Fig. [2Jd shows the Fermi 
surface that results from H adsorption on the understo- 
ichiometric surface. The suppression of the e 9 ' complex 
is weaker in the absence of Na on the surface, but still 
quite sufficient to push the e 9 ' pockets below E^. 

Although the Fermi surfaces in the stoichiometric unit 
cell are too complex to provide useful information, the 
position of the e s ' surface complex can nonetheless be 
gauged by examining the position of the e' g band center 
in comparison to the bulk. Note in Fig. [2^,, that the part 
of the e g ' band complex that gives rise to Fermi surface 
hole pockets (along T-M) decreases in energy more than 
the band complex as a whole, so a strong shift in band 
center will clearly indicate suppression of these pockets. 
We find for the under-stoichiometric, stoichiometric and 
over-stoichiometric surface the band shifts are -42, -120 
and -137 meV versus the bulk (two C0O2 layers beneath 
the surface) . We have also checked that the bulk numbers 
are unaffected by H adsorption on the surface (or by the 
existence of the surface itself) by comparing third-layer 
band centers to full bulk calculation band centers. Since 
the understoichiometric case exhibits the smallest band 
center shift and already has no surface e g ' hole pock- 
ets, these numbers demonstrate that the hole pockets 
are surely absent in the stoichiometric case. 
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Conclusions Formation of a surface results in the e g 
hole pockets shrinking and possibly disappearing. This 
effect is highly dependent on Na ordering and content 
on the surface. Surface hydroxyl termination, which can 
occur through hydrogen contamination in UHV, results 
in complete suppression of the e s ' hole pockets. Our 
energetics indicate that once hydrogen is adsorbed it will 
be very difficult to remove. This is in good agreement 
with hydrogen adsorption studies on other oxide surfaces. 
The Coulomb field of Na + or H + lowers the energy of 
the O electron bands which subsequently decreases Co- 
O bonding. Co-0 hybridization increases along the e g 
symmetry parts of the band complex that result in hole 
pocket formation, relative to the a\ g band complex. This 
results in e g hole pocket suppression. The decrease in 
size and absence of surface e g ' hole pockets is due to a 
combination of surface effects and surface contamination 
which agrees well with ARPES data. 
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sociatcship Program. 



FIG. 2: Band structure for (a) H adsorbed on overstoichio- 
metric surface and Fermi surface of (b) H adsorbed on un- 
derstoichiometric surface for Na 1/ / 3 Co02. Large and small 



circles indicate Na intercalant 
respectively. 



character and e 9 ' character, 



The dramatic effects of H on the surface owe to two 
underlying physical effects. The first is a simple filling of 
the Co-derived band complex due to electron donation 
from the adsorbed H. This is, in reality, a very small ef- 
fect, but since the e g pockets are already at the very top 
of the e g ' band complex and are the first part of the spec- 
trum to fill upon electron addition, at higher Na content, 
even small amounts of H adsorbed on the surface could 
destroy the small Fermi surface pockets. The second and 
dominant effect involves the positive Coulomb field cre- 
ated by ionized H just above the negative O ions which 
serves to exaggerate the main surface effect: O bands are 
strongly lowered, which decreases e g antibonding and 
destroys the pockets. This effect is operative even if the 
adsorbate is H + rather than atomic H, as may be the case 
for dissociated water molecules in the vacuum. Because 
O anti-bonding with e g and &\ g states is not uniform, 
there are also noticeable changes in the dispersion along 
r-M (See Figs. Hi and b). 
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